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I.  INTRODUCTION 


Many  factors  enter  into  the  design  of  a  rocket  motor 
(required  thrust,  specific  impulse,  operation  envelope, 
susceptibility,  safety,  cost,  flight  time,  weight,  size) . 

It  is  the  job  of  the  engineer  to  take  the  requirements  and 
to  choose  a  propellant  and  motor  configuration  that  best 
meet  these  needs.  In  order  to  do  this,  predictions  or 
calculation  of  performance  must  be  made.  Because  of  the 
high  cost  of  test  firing,  it  is  desirable  to  use  computer 
codes  that  will  accurately  predict  the  theoretical 
performance  of  various  compounds,  thereby  eliminating  the 
need  to  conduct  test  firings  on  all  but  a  few  acceptable 
candidates.  In  recent  years  these  codes  have  become 
available  in  PC  versions. 

After  the  propellants  are  selected,  performance  losses 
must  be  estimated.  With  these  losses,  actual  pressure- time 
and  thrust -time  profiles  and  specific  inpulse  can  be 
calculated  for  various  solid  propellant  grain  designs  and  or 
ramjet  configurations. 

It  was  the  intent  of  this  investigation  to  modify  two 
existing  computer  programs  to  cover  a  wider  variety  of 
performance  and  design  problems.  Specifically,  MICROPEP 
[Ref.  1]  (an  equilibrium  thermochemical  "propellant 
evaluation  code"  developed  by  the  Naval  Air  Warfare  Center 
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Weapons  Division,  China  Lake,  California)  was  to  be  modified 
to  allow  input  of  the  desired  nozzle  area  ratio  and  ambient 
pressure  so  that  effects  of  non- ideal  expansion  could  be 
accurately  calculated.  The  performance  value  for  shifting 
equilibrium  flow  to  the  nozzle  throat  and  frozen  composition 
flow  afterwards  is  often  a  good  approximation  to  the  actual 
flow  and  was  also  to  be  incorporated. 

In  order  to  incorporate  airbreathing  applications  the 
inlet  air  temperature  and  composition,  motor  static 
pressure,  nozzle  contraction  ratio  and  nozzle  discharge 
coefficient  are  required  inputs.  With  these  inputs  the 
program  should  calculate  the  vitiated  air  composition  and 
heat  of  formation  as  well  as  the  actual  combustion  chamber 
Mach  number  and  stagnation  pressure  and  the  associated 
performance  values. 

ROCKET  is  a  solid  propellant  rocket  motor  internal 
ballistics  code  orginally  developed  by  the  Lockheed 
corporation.  With  igniter,  propellant  and  insulation 
geometry  and  physical  porperties  input  it  calculates  the 
pressure- time  and  thrust -time  profiles.  Some  of  the 
subroutines  have  never  worked  correctly  and  many  unused 
portions  of  the  code  had  never  been  removed.  Docximentation 
for  the  code  is  also  very  limited.  Modifications  to  ROCKET 
were  required  to  make  the  nozzle  loss  coefficient  an  input, 
to  update  the  basic  algorithm  so  that  multiple  grains  could 
be  burned  at  the  same  time  and  to  incorporate  plot  routines 
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for  Che  pressure- time  and  thrust -time  profiles.  Average 
thrust,  average  pressure  and  total  impulse  were  also  desired 
output  values. 
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II.  MICROPEP 


A.  BACKGROUND 

The  MICROPEP  program  is  a  PC  based  computer  program 
developed  by  the  Naval  Air  Warfare  Center  Weapons  Division, 
China  Lake,  California  and  is  used  to  estimate  the 
performance  of  solid  and  liquid  propellant  rockets.  It  also 
has  the  capability  of  determining  equilibrium  and  thermo¬ 
chemical  properties  at  a  specified  pressure- temperature 
point . 

MICROPEP  can  be  broken  into  three  parts.  The  first  part 
provides  an  evaluation  of  combustion  product  composition  and 
chemical  properties  assuming  constant  pressure,  equilibrium 
adiabatic  combustion.  The  unmodified  program  allows  input 
of  up  to  ten  ingredients  and  the  output  consists  of  the 
following  combustion  chamber  properties:  temperature, 
pressure,  enthalpy,  entropy,  ratio  of  specific  heats, 
chemical  composition,  molecular  weight,  heat  content  and 
speed  of  sound. 

The  second  part  of  the  program  provides  an  evaluation  of 
the  throat  and  exhaust  compositions  and  chemical  properties 
assuming  isentropic,  adiabatic  expansion  through  the  nozzle. 
The  unmodified  program  calculates  these  properties  for  two 
separate  conditions: 
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1.  shifting  equilibrivun  flow  -  assumes  that  the 
reaction  kinetics  are  fast  enough  to  maintain  chemical 
equilibrium  throughout  the  expansion  process 

2 .  frozen  composition  flow  -  assumes  that  the  reaction 
rates  are  so  slow  that  the  exhaust  composition  is  the  same 
as  the  chamber  composition.  The  output  for  the  throat  and 
exhaust  are  the  Scime  as  for  the  chamber. 

The  third  part  of  the  program  provides  an  evaluation 
of  nozzle  performance  parameters  assuming  one -dimensional 
flow,  Mach  1  at  the  throat  and  exit  pressure  equal  to 
cimbient  pressure.  Output  consists  of  the  following  for  both 
frozen  and  shifting  flow:  specific  impulse,  ratio  of 
specific  heats  for  the  chamber  to  throat  process  (the  so- 
called  "process  gamma"),  characteristic  exhaust  velocity, 
nozzle  expansion  ratio,  density  specific  impulse,  vacuum 
impulse,  ratio  of  nozzle  throat  area  to  mass  flow  rate  and 
nozzle  thrust  coefficient.  An  example  of  the  output  from 
the  code  is  shown  in  Table  1.  It  should  be  noted  that  the 
units  are  a  mix  between  English  and  SI. 

B.  MODIFICATIONS 

None  of  the  modifications  made  to  the  program  used  any 
new  solution  optimization  algorithms.  The  existing 
subroutines  were  used,  as  needed,  to  affect  the  changes.  A 
flow  diagram  of  the  MICROPEP  program  is  shown  in  Figurel, 
with  the  modifications  shown  in  bold.  The  changes  made  are 


as  follows: 


Table  Is 

NEWPEP  SAMPLE  OUTPUT 


1  ****  NEWPFP  -  Fob.  1990  *♦** 

♦  XX  *  05/16/93  *  DU  ♦*  OEMS  ****  COMPOSITION  ******* 

OXYGEN  (GAS)  0  0.00001  20 

HYDROGEN  (GASEOUS)  0  0.00001  2)) 

INGREDIENT  WEIGHTS  (IN  ORDER)  AND  TOTAI.  WEIGHT  (T.AST  ITEM  IN  I.IST) 

94.0000  6.0000  100.0000 

THE  PROPELLANT  DENSITY  IS  0.0000.1  LB/CU-IN  OR  0.0003  GM/CC 
NUMBER  OF  GRAM  ATOMS  OF  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 
,  5.952381  H  5.875000  O 
****t**»t****t***n*nttttt  ‘CHAMBER  RE.SULTS  FOr.T,OW‘ 

T(K)  T(F)  P(A!t4)  P(PSI)  ENTHAl.PY  ENTROPY  CP/CV  SGAMtlA  RT/V 

3272.  5430.  34.01  500.00  0.00  293.21  1.2069  1.1361  7.204 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  4061.859  AND  4061.860  FT/SEC 

SPECIFIC  HEAT  (f40IjAR)  OF  GAS  AND  TOTAL-  11.590  11.590 

NUMBER  MOLS  GAS  AND  CONDENSED-  4.6693  0.0000 

2.63700  H20  1.32936  02  0.45880  HO  0.11383  O 

8.93E-02  H2  3 . 77E  02  H  3.32E-03  H02  6.49E-06  03 

THE  MOLECULAR  WRIGHT  OF  THE  MIXTURE  IS  21.417 

TOTAL  HEAT  CONTENT  (298  REF)  -1722.111  CAL/GM 
SENSIBLE  HEAT  CONTENT  (298  REF) -1 393 . 949  CAT./GM 

»  «  *  *  »  “  4  •  t  “  *  <  4  I  •  «  »  *  >  «  t  <  »  «  t  I  exhaust  results  F0LI.0W“  •  4  4  4  4  4  4  4  ...  4  I  4  4  I  4  4  4  4  4  4  4  4 

T(K)  T(F)  P(ATTl)  P(rsi)  ENTHAl.PY  ENTROPY  CP/CV  SGAMMA  RT/V 

2192.  3486.  1.00  14.70  -86.91  293.21  1.2123  1.1733  0.225 

DAMPED  AND  UNDAMPED  .SPEED  OF  SOUND-  3252.018  AMD  3252.018  FT/SEC 

SPECIFIC  HEAT  (MOLAR)  OF  GAS  AMD  TOTAL-  11.349  11.349 

NtR'tBER  MOLS  GAS  AND  CONDENSED-  4.4479  0.0000 

2.93879  H20  1.43343  02  O.0G25S  HO  0.00648  o 

5.62E-03  H2  P.50E-n4  H  1.62E-04  H02 

THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  IS  22.48.3 
TOTAL  HEAT  CONTENT  (298  REE)  -1131.643  CAL/GM 
SENSIBLE  HEAT  CONTENT  (298  REE)-  819.350  CAL/GM 
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Table  1:  (cont.) 

MEWPEP  SAMPLE  OUTPUT 


****** ***♦pERFOR^tANCE:  FROZEN  ON  FIRST  LINE,  SHIFTING  ON  SECOND  LINE******* 

All  exact  method  for  determining  throat  conditions  was  used 
The  frozen  &  shifting  STATE  gammas  for  the  throat  are:  1.2063  1.1360 

ISentroplc  Exponent  shown  below  is  the  gamma  for  the  chamber  to  throat 
PROCESS . 

lMPUf.SE  IS  EX  T*  P*  C*  ISP*  OPT  EX  D  ISP  A*M.  EX  T  ADI! 

263.3  1.2089  2981.  19.83  5691.9  5.09  0.1  0.35390  1728.  581778. 

275.0  1.1360  3096.  19.63  5815.9  222.9  5.77  0.1  0.36161  2192.  770299. 

**************************** *throat  results  FOI.LOW* ************************* 

T(K)  T(F)  P(ATM)  F(PS1)  ENTHALPY  ENTROPY  CP/CV  SGAMHA  RT/V 

3096,.  5114.  19.63  288.61  -16.14  293.21  1.2063  1.1360  4.250 

DAMPED  AND  ONDTVMPED  SPEED  OF  SOHNO-  3929.313  AND  3929.314  FT/SEC 

SPECIFIC  HEAT  (MOLAR)  OF  GAS  AND  TOTAL’  11.621  11.621 

NUMBER  MOLS  GAS  AND  CONDENSED’  4.6196  0.0000 

2.70089  H20  1.35026  02  0.38104  HO  0.08796  O 

7.01E-02  H2  2.71E-02  H  2.28E-03  H02  3.74F-06  03 

THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  IS  21.647 

TOTAL  HEAT  CONTENT  (298  REF)  ’1618.754  CAL/GM 
SENSIBLE  HEAT  CONTENT  (298  REF) -1297 . 611  CAL/GM 

****  NEWPEP  -  Feb.  1990  **** 
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Figure  1  : 

PEP93  Flow  Diagram 


PEP93  Flow  Diagram 
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Call  SBAL 

Calculate  exhaust  conditions 
for  equilibriunn  flow  and 
contant  entropy 


Guess  new  exit 
pressure 

T 


Call  DESIGN 

Calculate  throat  conditions 
assuming  Mach  1 


Figure  1 ;  (cont.) 
PEP93  Flow  Diagram 
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1.  The  maximum  allowable  number  of  input  ingredients 
has  been  increased  to  fifteen  and  the  gas  constant  for  the 
mixture  is  calculated. 

R  =  R/M 

2.  The  original  program  was  extremely  concise,  leaving 
off  units.  Choosing  clarity  over  brevity,  units  were  added 
for  all  data. 

3.  Performance  calculations: 

a.  Shifting  equilibrium  versus  frozen  composition 
flow-  The  theoretical  performance  of  a  rocket  engine  varies 
with  the  two  assumptions.  If  the  reaction  rates  of  the 
chemicals  are  fast  enough  to  maintain  equilibrium  throughout 
the  nozzle  then  the  additional  energy  gained  by  the 
reassociation  of  molecules  and  atoms  will  be  reflected  in  an 
increased  kinetic  energy  at  the  exhaust  nozzle  exit. 

Looking  at  the  thrust  and  specific  impulse  equations 
F  =  mVg/gc  -  (Pe-Po)Ae 
Isp  = 

it  is  apparent  that  theoretical  performance  would  be  highest 
for  shifting  flow  and  lowest  for  frozen  flow.  Since  actual 
performance  has  been  measured  between  these  two  conditions, 
modifications  have  been  made  to  MICROPEP  to  calculate 
shifting  flow  to  the  throat  and  frozen  flow  afterwards. 

This  accounts  for  the  higher  reaction  rates  due  to  higher 
temperatures  and  pressures  coupled  with  a  flow  Mach  number 
less  than  one  prior  to  the  throat.  After  passing  through 
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the  throat  the  flow  accelerates  past  Mach  l  and  the  pressure 
drops  rapidly  to  near  ambient  pressure,  making  it  unlikely 
that  the  reaction  kinetics  will  keep  up  with  the  flow 
velocity.  Therefore,  the  assumption  of  frozen  flow  is 
reasonable. 

b.  Non- Ideal  Expansion  -  It  can  be  easily  shown 
that  the  maximum  thrust  occurs  when  the  exit  pressure  is 
equal  to  ambient  pressure  (Ref.  2] 

F  =  mVg/g^  -  (Pe-Po)Ae 

Taking  the  derivative  of  thrust  with  mass  flow  rate 
and  ambient  pressure  constant  gives 

dF  =  mdVg/g^  +  (p^-p^)6A^  +  A^dpg 
Using  the  continuity  equation  and  the  momentum  equation  for 
ideal  flow 

dF  =  (pAV)gdVe  +  (Pg-Po)dAe  +  A^(-pVdV)^ 
dF  =  (Pe-Po)dAe 

therefore 

dF/dAg  =  0  if  Pe=Po 

In  an  actual  rocket  motor  the  exit  area  is  constant,  so  the 
expansion  will  only  be  ideal  for  one  operating  altitude.  An 
option  was  added  to  input  nozzle  expansion  ratio  and  the 
ambient  pressure.  Performance  values  for  non- ideal  expansion 
can  now  be  calculated  for  different  ambient  pressures. 
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4 .  Exhaust  kinetic  energy 

One  of  the  output  values  of  NEWPEP  was  labeled 
adiabatic  head.  This  value  is  calculated  using  the  formula 


and  is  therefore  labeled  kinetic  energy  of  exhaust. 

5.  Total  and  sensible  heat  content 

These  were  included  in  the  original  NEWPEP  but  the 
definitions  were  not  clear.  The  present  values  are: 

o  o 

a.  Total  Heat  Content  =  Epn£hfi  -  Epn^^h^^ 

(without  water  condensation) 

o  o 

b.  Sensible  Heat  Content  =  Epn£hfi  -  Epn^h^^ 

(with  water  condensed) 

6.  Air-breathing  motors 

a.  Vitiated  air 

The  theoretical  performance  of  the  combustor 
varies  with  the  composition  and  temperature  of  the  incoming 
air.  The  vitiated  air  option  requires  the  user  to  select 
the  vitiator  ingredients,  inlet  stagnation  temperature  and 
pressure.  Using  this  information  a  pressure- temperature 
point  calculation  is  made  for  the  inlet  composition  and  heat 
of  formation.  This  "vitiated"  air  is  chen  automatically 
entered  as  a  user  defined  ingredient  for  the  combustor 
calculations . 
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b.  Combustor  Mach  number 

As  combustor  Mach  number  increases  the 
difference  between  static  and  stagnation  pressure  increases 
Pt/p  =  (l+(  (7-1)/2)*m2)"(y/y-1) 

The  unmodified  program  assumes  the  Mach  number  is  zero  and, 
therefore,  stagnation  pressure  is  equal  to  static  pressure. 

To  account  for  finite  Mach  numbers,  the  nozzle 
contraction  ratio  can  be  input.  If  this  option  is  chosen 
the  program  first  calculates  the  nozzle  and  exhaust 
conditions  for  a  chamber  Mach  number  equal  to  zero.  Then, 
using  the  calculated  process  gamma  and  given  nozzle 
contraction  ratio,  the  combustion  chamber  Mach  number  is 
solved  for  using  an  interval  halving  method.  Using  a  lower 
bound  of  M  =  0  and  an  upper  bound  of  M  *  1.001,  the 
following  equation  is  iterated  until  convergence  is  reached. 

0  =  M*{A4/A5)  -  [2/(7+l)  +  (7-l)/(7+l)*M2)]^t  (7  +  1)/2{7-1)] 
The  bisection  method  was  chosen  because  of  its  guaranteed 
convergence  and  the  small  range  of  Mach  numbers  to  be 
considered. 

c.  Nozzle  discharge  coefficient 

One  of  the  basic  assumptions  of  MICROPEP  is  one- 
dimensional,  isentropic  flow.  In  order  to  more  accurately 
predict  mass  flow  rate  through  the  nozzle,  a  streamline 
curvature  correction  (C^j)  is  now  an  allowable  input.  With 
the  nozzle  discharge  coefficient  the  continuity  equation 
becomes 
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mg  =  PtAtCd[7/RTt]°-^[2/(7+l)]"t(7+l)/2(7-l)3 
=  Pt^tCd/C* 

If  an  accurate  value  is  not  known  the  program  default  is 

1.0. 

7.  Equivalence  Ratio 

One  definition  of  equivalence  ratio  is  the  ratio  of 
oxidizer  required  for  complete  combustion  of  the  fuel  (to 
the  highest  oxidation  state)  to  that  available  for 
combustion  [Ref.  3] .  Using  this  definition,  subroutine 
EQUIVALENCE  solves  the  following  equation 

4>  =  ni  X  (no)c]/no 

where  K  is  niimber  of  reactant  elements, 

n^  is  the  number  of  moles  of  fuel  element  i, 
in^)  ^  is  the  number  of  moles  of  oxygen  required  for 
complete  combustion  of  one  mole  of  fuel  element  i, 

n^  is  the  number  of  moles  of  oxygen  available. 

Note:  The  program  attempts  to  account  for  all  common  fuel 

species.  See  Appendix  C:  Subroutine  EQUIVALENCE  for  a 
complete  list  of  the  elements  incorporated. 

C.  REVISED  INPUT  AND  OUTPUT 

Appendix  A  contains  the  revised  "users  guide"  for  PEP93. 
A  sample  output  is  shown  in  Table  2 . 
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Table  2  s 

PEP 9 3  SAMPLE  OUTPUT 


****  NEWPEP  - 

Jan. 

1993  *♦** 

*  XX  * 

INGREDIENT 

MASS 

HF 

DENSITY 

COMPOSITION 

(gm) 

(kcal/kg) 

(kg/m**3) 

OXYGEN  (GAS) 

94.00 

.0 

.0000 

20 

HYDROGEN  (GASEOUS) 

6.00 

.0 

.  0000 

2H 

VOLUME  PERCENT  OF  INGREDIENTS  (TN  ORDER) 

94.0000  6.0000 

THE  PROPELLANT  DENSITY  IS  .00001  LB/CU-IN  OR  .0003  GM/CC 
THE  EQUIVALENCE  RATIO  IS  ■ 5066 

NUMBER  OF  GRAM  ATOMS  OF  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 
5.9523B1  H  5. 076000  O 


*♦*******♦«**♦♦♦*********  ‘  * ‘CHAMBER  RESm.TS  FOt,I,OW 


TEMP  PRESSURE 

(K)  (MPa/ATM/PSl) 

3271.9  3.447/  34.02/  500.00 


EUTHAf.Py 

(kJ/kg) 

.0000 


ENTROPY  CP/CV 
(kJ/kg-K) 

12.260  1.2069 


SGAMMA  P.i/itJ 
(MPa/kmol ) 


1.136)  73031.150 


DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  1230.061  AND  1230.062  m/ppc 

SPECIFIC  HEAT  (MOLAR)  OF  GAS  7VNO  TOTAL-  404.92  484. 122  J/):g  K 

NUMBER  MOIiS  GAS  AND  CONDENSED^^  4.6693  .0000 


(♦-liquid,  li-gol  id) 

2.63700  1120  1.32935  02 

0.93E-02  112  3.77F  02  H 

THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE 
THE  GAS  CONSTANT  IS  308.22  J/kq  K 

TOTAL  HEAT  CONTENT  (290  REF) 
SENSIBLE  HEAT  CONTENT  (298  REF) - 

4. EXHAU.ET 


.45081  HO  .11383  O 

3.32R-03  H02  6.49E06  03 

IS  21.417  qm/niolp 

7205.201  kJ/kg 
5032.25?  k,I/);g 

RESULTS  FOLLOW . *.♦♦*»♦««♦**♦* 


TEMP 

(R) 

2)91.9 


PRESSURE 
(MPa /ATM/ PS I) 
.101/  1.00/  14.70 


ENTHAl-Py 

(kJ/kg) 

.3636.2370 


ENiRory  CP/CV 
(kJ/):g-K) 

12.268  1.2123 


SGAI-1MA  Pl/tii 
(MPa/ktnol ) 


1.1733  2278.684 


DAMPED  AND  UNDAMPED  SPEED  OF  SOUND^  991.215  AIHl  991.215  tn/sec 

SPECIFIC  HEAT  (MOLAR)  OF  GAS  AND  TOTAL-  474.02  474.822  J/kg  K 

NUMBER  MOLS  GAS  AND  CONDENSED-  4.4479  .0000 
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Table  2s  (cont.) 
PEP93  SAMPLE  OUTPUT 


(*=llguid,  5t=solld) 

2.93879  1120  1.43343  02  .06255  HO  .00647  0 

5.62E-03  112  8.50E-04  H  1.62E-04  H02 

THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  IS  22.483  gm/mole 
THE  GAS  CONSTAMT  IS  369.81  J/kg-K 

TOTAL  !IEAT  CONTEIIT  (299  REF)  -  4734.795  kj/kg 

SENSIBLE  HEAT  CONTENT  (298  REF)-  3428.162  kJ/kg 

ttt***t*******t***********»t********»********t**»***********************t***t 

All  exact  method  for  determining  throat  conditions  was  used 

Tlie  frozen  &  shifting  STATE  gammas  for  the  throat  are:  1.2063  1.1360 

GAtWA  NU  shown  below  Is  the  gamma  for  tlie  chamber  to  t)iroat  PROCESS. 

»*********PErFORMANCE:  FROZEN  ON  FIRST  LINE,  SHIFTINO  ON  SECOND  LINE, <***** 

*♦♦***«***. 4**gnTFriNG  TO  THROAT/FROZEN  AFTERWARDS  ON  THIRD  I.TNE* ********** ♦ 


SPECIFIC  GAMMA 
IMPULSE  NU 

T» 

P* 

C* 

TSP* 

Ae/A* 

D  ISP 
(gm- s/ 

A*/m 

(cm**2/ 

Te  Cf 

(sec) 

(K) 

(MPa) 

(m/s) 

(sec) 

cm* *3) 

kg/s) 

(K) 

263.3  1.2000 

2962. 

1.930 

1733.3 

5.099 

2.017 

5.0272 

1728.1.4895 

275.0  1.1360 

3096. 

1.990 

1772.7 

222.9 

5.768 

2.107 

5.1415 

2192.1 .5215 

267.2  1,1360 

3096. 

1.990 

1772.7 

222.9 

5.110 

2.047 

5.1115 

1016.1.1781 

**^***nt***********tt*»*t**  *1HR0AT  RR.SULTS  FOLLOW*  «♦♦«*♦*♦*****«♦***♦♦♦**♦♦ 

TEMP  PRESSURE  ENTHALPY  ENTROPY  CP/CV  .SGAMt-lA  Pj/ni 

(K)  (MPa/ATM/PSl)  (kJ/kg)  (kJ/kg-K)  (MPa/kmol) 

3096.5  1.990/  19.64/  288.61  -675.5011  12.268  1.2334  1.1441  2193.979 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  1211.075  AND  1211.075  m/s°r 

SPECIFIC  HEAT  (MOLAR)  OF  GAS  AND  TOTAL--  4  39.30  4  39.299  J/fcg  K 

NU14BER  MODS  GAS  AND  CONDENSED-  4.6196  .0000 

( *  =1  Iquld,  Si-sol  id) 

2.70090  (120  1.35026  02  .38103  HO  .08795  O 

7.01B-02  H2  2.71E  02  H  2.28E  03  H02  3.74E  06  03 

THE  MOLECUIjAR  WEIGHT  OF  THE  HIXIUPE  IS  21.647  gm/mole 
THE  GAS  CONSTANT  IS  384.09  J/kg-K 

TOTAL  HEAT  COiriEHT  (298  REF)  -  6772.908  kJ/kg 

SENSIBLE  HEAT  CONTENT  (290  REF)-  5429.259  kJ/kg 

FROZEN  6  SHIFTING  KINETIC  ENERGY  OF  EXHAU.ST  504773.  770295.  m**2/s**2 
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III.  ROCKET 


K.  BACKGROUND 

1.  ROCKET  is  a  PC  based  program  originally  written  by 
the  Lockheed  Corporation  and  then  modified  by  the  Naval  Air 
Warfare  Center  Weapons  Division.  The  intent  of  the  original 
program  was  to  calculate  the  pressure- time  and  thrust -time 
profiles  of  a  multiple  grain  rocket  motor  using  the 
following  assumptions: 

a.  Combustion  products  are  ideal  gases. 

b.  Burning  rate  follows  Veille's  Law 
r  =  ap^” 

c.  Effects  of  mass  addition  and  erosive  burning 
are  negligible. 

d.  Inertia  of  the  chamber  gases  is  negligible. 

e.  Variation  of  C*  with  pressure  is  described  by 

C*  =  [Pref/lOOO]X 

f.  Propellant  burning  rate  temperature 
sensitivity  is  described  by 

r  =  f^Qpexp  [0  . OlcTp  (T^- 70)  ]  ,  where  ffp  =  (l-n)I^ 

g.  Properties  of  the  chamber  gases  are  mass 
averages  of  the  individual  gases. 

2.  Allowable  inputs  were: 

a.  Number  and  type  of  grains 
(1)  Igniter 
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(2)  Insulation 

(3)  Propellant 

(4)  Igniter  Motor  (for  large  rockets) 
b.  Motor  parameters 


(1) 

Nozzle 

throat  and  exit  areas 

(2) 

Total  : 

motor  volume 

(3) 

Ambient  temperature  and  pressure 

(4) 

Throat 

closure  blowout  pressure 

(5) 

Initial  motor  pressure 

(6) 

Nozzle 

discharge  coefficient  and  half 

angle 

correction 

(7) 

Ratio 

of  specific  heats 

(8) 

Throat 

area  design  pressure 

(9) 

Throat 

erosion  rate (radial) 

c.  Grain  Parameters 

(1)  Burning  rate 

(2)  Burning  rate  exponent 

(3)  Temperature  sensitivity 

(4)  Reference  pressure  for  C*  correction 

(5)  Characteristic  exhaust  velocity 

(6)  Pressure  correction  exponent  for  C* 

(7)  Density 

(8)  Ignition  time 

(9)  Flame  spread  ignition  delay  time 

(10)  Burning  surface  area 
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3.  Using  these  inputs  the  program  calculates  chamber 
pressure  by  integrating  the  following  mass  conservation 
equation: 

dpc/dt  =  (I/Wq)  [RT{rppAbapc"-p^h9c/C*) 3 
where  is  the  instantaneous  gas  volume.  Then  the  thrust 
coefficient  is  calculated  using 
Cp  =  {272/(7-1)*(2/(7+1)  )"[(7+l)/(7-l)] 

*  [1- (Pg/Pc)"  (7-1)/7]  +  (Pe'Pamb^Ae/ (PcCd^th) 

With  the  values  of  pressure  and  thrust  coefficient  known, 
thrust  was  solved  for  using 
F  =  CpP^h>^Cd 

B.  MODIFICATIONS 

A  flow  diagram  of  the  program  is  shown  in  Figure  2  and  a 
modified  user's  guide  is  given  in  Appendix  B. 

1,  Calculation  of  actual  thrust  and  specific  impulse. 
Several  factors  contribute  to  the  reduction  in  theoretical 
thrust  and  specific  impulse  from  the  theoretical  values 
(F' ,  Isp' ) .  In  general  these  can  be  put  into  two  separate 
categories;  losses  due  to  incomplete  combustion  in  the 
chamber  and  nozzle  losses  (two-phase  flow,  divergence, 
etc.) .  The  ROCKET  program  has  been  modified  to  allow  these 
losses  to  be  accounted  for  using  two  efficiencies,  tjq*  and 
Vcf 

a.  Combustion  Loss  Coefficient 

tjq*  is  primarily  a  function  of  the  gas  residence 
time  and  is  determined  empirically  [Ref.  4] .  Typical  values 
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MAIN 

Controlling  program 

I 

Call  INPUT 

Read  in  general  motor 
description,  grain  parameters 


Calculate  propellant 
volume  and  weight 


Correct  C*for 
actual  chamber 
pressure 


NoJ 


Calculate 
throat  area 
based  on 
input 
pressure 


Input  throat  area 
assumed  correct  and 
program  only  runs  once 

S 

Call  EVES 

Integrate  pressure  equation  over 
small  time  increments 


Call  DIFEQ 

Calculate  throat  erosion,  change  in 
pressure,  new  chamber  pressure 
and  mass  flow  rate 


No 


Call  PRINT 

Print  current  pressure,  change  in 
pressure,  free  volume,  propellant  weight 
consumed,  insulation  weight  consumed 
and  throat  area 


Yes- 


Call  GETAVG 

Calculate  average  pressure 
over  total  bum  time 

_ 

> 

Print  average  pressure 
and  maximum  pressure 

1 

Call  THRUST 
Calculate  thrust, 
coefficient  of  thrust, 
specific  impulse  and  total 
impulse 


Call  GETAVG 
Calculate  average  thrust 
over  total  bum  time 


Print  average  thrust 
and  total  impulse 


Call  PLOT 

Plot  thrust  vs  time  and 
pressure  vs.  time 


END 


Figure  2: 

ROCKET  Flow  Diagram 
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for  aluminized  solid  propellants  are  between  .96  and  .995 
for  residence  times  greater  than  10  seconds.  For  these 
propellants  drops  off  rapidly  for  residence  times  below 
10  seconds. 

b.  Nozzle  Loss  Coefficient 

This  can  be  expressed  by  tRef .  5] . 

T/cf  =  1-  (ediv+^kin+^bl  +  ^sub'^^eros+^tp) 

where 

^div  =  half  angle  correction  for  conical  nozzles 
^kin  ~  correction  for  shifting  equilibrium  vs. 

frozen  composition  flow 
ejgj  =  boundary  layer  loss 

6gub  ®  correction  for  flow  around  submerged 
nozzles 

6gros~  losses  due  to  nozzle  throat  erosion 
6j.p  =  velocity  and  thermal  lag  of  condensed 
species  in  chamber 

For  use  in  ROCKET  the  input  value  of  C*  should  be  C* 
theoretical  times  .  Allowable  inputs  of  nozzle 

discharge  coefficient  and  half -angle  correction  have  been 
replaced  by  the  nozzle  loss  coefficient,  and  thrust  is 
calculated  using 

2.  Plotting  routines 

A  plotting  subroutine  has  been  added  which 
automatically  takes  the  time  and  thrust  values,  determines 
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the  maximum  value  of  thrust  and  plots  the  thrust -time 
profile  scaled  to  maximum  thrust.  The  same  is  done  for  the 
pressure- time  profile.  In  addition,  time,  pressure  and 
thrust  are  printed  to  a  separate  data  file. 

3.  Total  impulse  is  calculated  using  the  following 
formula: 

I-ji  =  Ftjj 

where 

F  is  the  average  thrust 
tj3  is  the  total  burn  time. 

4.  Improved  efficiency 

Since  the  original  program  was  written  by  Lockheed 
for  use  on  a  mainframe  computer,  ROCKET  has  undergone  some 
modifications.  These  changes  resulted  in  a  program  that  was 
inefficient  because  several  input  variables  were  undefined 
in  the  users  manuals.  In  order  to  increase  efficiency,  more 
than  25  variables  and  over  300  lines  of  code  associated  with 
these  variables  were  removed.  The  result  is  a  program  that 
still  gives  the  output  discussed  above  and  does  indeed  work 
well  for  more  than  one  grain. 

C.  REVISED  INPUT  AND  OUTPUT 

A  sample  input  file  is  given  in  Appendix  B  with  a  rocket 
motor  description  shown  in  Figure  3.  A  sample  output  file 
is  shown  in  Table  3. 


23 


propellant  grain  2 


Table  3i 

ROCKET  SAMPLE  OUTPUT 


IBpercent 

Motor  Performance  Program  IBM- PC  Version  1.0 
INPUT  DATA. . . 

Normal  Input  Is  finished 
General  Configuration  PARAMETERS  are  ... 


Throat  AREA 
Exit  AREA 
Expansion  RATIO 
Total  Motor  VOLUME 
AMBIENT  Temperature 
AMBIENT  Pressure 
Closure  BI^OWOUT 
Pzero 

CP  Efficiency 
Gamma 

Throat  DESIGN  Pressur 
Tliroat  EROSION 
MOTOR  HAS  4  GRAINS 
DESCRIPTION  OF  GRAIN  1 
BtJRH  RATE 
BURN  RATE  EXP. 

PI  StIB  K 

BtfRN  REF.  PRESS. 

C  STAR 
C  STAR  EXP. 

DENSITY 
IGNITION  TIME 
DELTA  IGN.  TIME 


8.6577 
105.4500 
12.1799 
8042.5 
70.0 
10.1 
35.0 
14.7 
.9606 
1 . 1443 
1030.0 

.0005 


.5800 

.3500 

.2000 

1030.0 

5106.7 

.0450 

.0623 

.0000 

.0000 


WEB 

.000 

.301 

.762 

1.143 

1.524 

1.905 

2.286 

2.666 

3.047 

3.428 

3.809 


BtJRN  AREA 

744 . 150 
759.750 

771.250 
778.650 
781.950 

781.150 

776.250 

767.250 
754.100 
736.900 
715 . 600 


Table  3:  (cont.) 
ROCKET  SAMPLE  OUTPUT 


DESCRIPTION  OF  GRAIN  2 

BURN  RATE 

.5800 

BURN  RATE  EXP. 

.3500 

n  SUB  K 

.2000 

BURN  REF.  PRESS 

1030.0 

C  STAR 

5106.7 

C  STAR  EXP. 

.0150 

DENSITY 

.0623 

IGNITION  TIME 

.0000 

DELTA  IGN.  TIME 

.0000 

WEB 

BURN  AREA 

.000 

744 . 150 

.381 

759 . 750 

.762 

771.250 

1.143 

778.650 

1.524 

781.950 

1.905 

781.150 

2.286 

776.250 

3.666 

767.250 

3.047 

754.100 

3.420 

736.900 

3.009 

715.600 

DESCRIPTION  OF  GRAIN  3 

BURN  RATE 

.5800 

BURN  RATE  EXP. 

.  3500 

n  SUB  K 

.2000 

BURN  REF.  PRE.SS 

1030.0 

C  STAR 

5106.7 

C  STAR  EXP. 

.0150 

DENSITY 

.0623 

IGNITION  TIME 

.  0000 

DETjTA  IGN.  TIME 

.0000 

THI.S  GRAIN  IS  INSULATION  WITH 

AN  ABLATION  RATE  OF 

WEB  1 

BURN  AREA 

.000 

1488.300 

.381 

1519 . 500 

.  762 

1542 . 500 

1 .143 

1557.300 

1.524 

1563.900 

1.905 

1562 . 300 

2.286 

1552.500 

2.666 

1534.500 

3.047 

1500.200 

3 .420 

1473 . 800 

3.809 

1431 .200 

.05000 
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Table  3s  (cont.) 
ROCKET  SAMPLE  OUTPUT 


DESCRIPTION  OF  GRAIN  4 
BURN  RATE 
BURN  RATE  EXP. 

PI  SUB  K 

BURN  REF.  PRE.SS. 

C  STAR 
C  STAR  EXP. 

DENSITY 
IGNITION  TIME 
DELTA  IGN.  TIME 
THIS  GRAIN  IS  AN  IGNITER 
TIME  WDOT 


.0000 

1  . 

4000 

.2500 

1 . 

9000 

TOTAL 

PROPELI.ANT 

WEIGNT- 

362 .4 

TIME 

PC 

Poor 

FREE  VOL 

.000 

14.70 

22173.60 

221 7 . 59 

.001 

27.75 

27381.44 

2217.75 

.001 

38.26 

27371.20 

2217.88 

.001 

51.96 

29634.08 

2218.07 

.002 

66.64 

31490.92 

2218.27 

.002 

76.09 

32543.34 

2210.41 

.002 

87.46 

33467.93 

2218.55 

.003 

98.30 

34278.57 

2218 . 70 

.003 

109.38 

34986.84 

2218.86 

.003 

120.68 

35602 . 51 

2219.02 

.004 

132.16 

36134.02 

2219 . 18 

.004 

143.80 

36588.74 

2219 . 35 

.004 

155.57 

36973.17 

2219 . 52 

.005 

167.45 

37293.06 

2219 . 70 

.005 

191.40 

37759.36 

2220.06 

.006 

215.74 

38023 .35 

2220.43 

.007 

240.11 

38114.57 

2220.02 

.007 

264.49 

38057.92 

2221.22 

.008 

288.80 

37874.59 

2221.63 

.009 

312.95 

37582.82 

2222.05 

.009 

336.89 

37198.40 

2222.48 

.010 

360.55 

36735.09 

2222.92 

.010 

383.89 

36204.92 

2223 . 37 

.011 

406.88 

35618.45 

2223 . 83 

.012 

451.65 

34312.61 

2224.77 

.014 

494.66 

32879.32 

2225.74 

.015 

535.78 

31366.61 

2226.73 

.016 

574.94 

29811.89 

2227.74 

.018 

612.09 

28244 . 19 

2228.78 

.019 

647.24 

26605.84 

2229.83 

.  020 

680.42 

25153.92 

2230.91 

.  021 

711 .65 

2.4661 . 14 

2232.00 

1.0000 

-0000 

-0000 

1000.0 

5106.7 

.0450 

.0623 

.0000 

.0000 


IGNITER  WEIGIITr- 


P  WEIGHT 

WDOTP 

.00 

13.55 

.01 

16.58 

.01 

18 . 39 

.02 

20.31 

.03 

22.03 

.04 

23,09 

.05 

24.09 

.06 

25.04 

.06 

25.94 

.07 

26.80 

.08 

27.62 

.09 

28.41 

.10 

29 . 16 

.11 

29 . 89 

.13 

31.26 

.15 

32.53 

.17 

33.72 

.19 

34.83 

.22 

35.88 

.24 

36,86 

.26 

37 . 79 

.29 

38.67 

.31 

39.50 

.34 

40.28 

.39 

41  .  73 

.41 

43.01 

.50 

44.22 

.56 

45.30 

.62 

4  6.27 

.69 

47.16 

.74 

47.97 

.90 

49.71 

IME  WT 
.00 
.00 
.00 
.01 
.01 
.01 
.01 
.01 
.01 
.02 
.02 
.02 
.  02 
.02 
.02 
.03 
.03 
.03 
.04 
.04 
.04 
.05 
.05 
.05 
.06 
.06 
.07 
.00 
.00 
.09 
.09 
.10 


.41 

MT>Or  AT 

.00  8.65770 
.00  8.55770 
2.42  8.65770 

3.24  8.65770 
4.11  0.65771 
4.71  8.65771 
5.32  8.65771 
5.95  8.65771 
6.59  8.65771 

7.24  8.65772 

7.90  8.65772 
8  56  8.G5V72 
9.23  8.65772 

9.90  0.65772 
11.25  8.65773 
12.61  8.65773 
13.97  8.65773 
15.32  8.65771 
16.66  8.65774 
17.99  8.65774 
19.30  8.65774 
20.59  8.65775 
21.86  8.65775 
23.11  8.65775 
25.53  8.65776 
27.85  8.65777 
30.06  8.65777 
32.15  8.65778 
34.13  8.65779 
36.00  8.65779 
37.76  0.65700 
39.42  8.65781 
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Table  3i  (cont.) 

ROCKET  SAMPLE  OUTPUT 


TIME 

FC 

FOOT 

FREE  VOL 

F  WRIGHT 

WDOTP 

INS  WT 

WDOT 

AT 

.023 

741.01 

22216.91 

2233.10 

.86 

49.39 

.  11 

40.97 

8.65701 

.024 

760.55 

20827.85 

2234.22 

.92 

50.01 

.11 

42.42 

8.65782 

.025 

794.35 

19498.47 

2235.34 

.99 

50.58 

.12 

43.78 

8.65783 

.026 

818.49 

16231.46 

2236.48 

1.05 

51 . 10 

.12 

45.05 

8.65783 

.029 

862 . 10 

15889 . 49 

2238.79 

1.19 

52.02 

.13 

4  7.34 

8.65785 

.032 

900.05 

13000.72 

2241 . 14 

1.32 

52.80 

.  15 

49.33 

8.65786 

.034 

932.96 

11954.15 

2243.51 

1.46 

53 .46 

.16 

51 . 05 

8.65707 

.037 

961.43 

10332.46 

2245.91 

1.59 

54.02 

.17 

52.54 

8.65789 

.039 

986.03 

0916.11 

2248 . 33 

1.73 

54.49 

.18 

53.82 

0.65790 

.042 

1007.24 

7684 . 14 

2250.77 

1.87 

54 .90 

.  19 

54.93 

8.65791 

.044 

1025.51 

6616.43 

2253.22 

2.01 

55.25 

.21 

55.88 

8.65793 

.047 

1041.23 

5693.39 

2255.69 

2.16 

55 . 54 

.22 

56.69 

8.65794 

.050 

1054.76 

4897.45 

2258 . 17 

2.30 

55.80 

.23 

57.40 

8.65795 

.052 

1066.40 

4212.11 

2260.66 

2.44 

56.02 

.24 

58.00 

8.65797 

.055 

1076.41 

3623.16 

2263.16 

2.59 

56.20 

.25 

58.52 

8.65798 

.060 

1092.43 

2683.66 

2260 . 17 

2 . 07 

56.51 

.28 

59 . 35 

8 . 65801 

.065 

1104.32 

1993.59 

2273.20 

3 . 16 

56.74 

.30 

59 . 97 

8.65803 

.073 

1116.71 

1280.38 

2280 . 79 

3.60 

56.98 

.  34 

60-61 

8 . 65807 

.003 

1126.01 

740.77 

2290.94 

4.19 

57 . 20 

.  38 

61 . 14 

8 . 65813 

.106 

1137.17 

266.38 

2313 . 06 

5.51 

57.46 

.49 

61.68 

8.65825 

.  183 

1148.30 

113.77 

2390.75 

9.94 

57.95 

.85 

62.25 

8.65865 

.249 

1155.79 

-2965.13 

2457.94 

13.81 

56.44 

1 .16 

62 . 64 

8.65900 

.255 

1140.96 

-2095.53 

2463.65 

14.14 

56 . 19 

1.19 

61.88 

8.65903 

.263 

1120.01 

-1334.14 

2471 .14 

14.57 

55.90 

1.22 

61  21 

0.65907 

.273 

1117.77 

-724.26 

2481.10 

15.15 

55.82 

1.27 

60.68 

8.65912 

.311 

1107.56 

-36.76 

2518.33 

17.29 

55.71 

1.45 

60.15 

8.65932 

.322 

1107.41 

2.71 

2528.26 

17.86 

55.73 

1.50 

60.14 

0.65937 

.532 

1110.03 

55.27 

2732 . 80 

29.61 

56.30 

2.49 

60.70 

8.66047 

.736 

1128.33 

42.24 

2934.31 

41.20 

56.01 

3.45 

61.24 

8.66154 

.982 

1138.41 

41.30 

3178.18 

55.22 

57.32 

4.62 

61.77 

8.66282 

1.228 

1148.53 

41.29 

3424.17 

69.37 

57.83 

5.80 

62.31 

0.66410 

1.576 

1158.84 

26.43 

3775.79 

89.60 

58.34 

7.48 

62.85 

0.66592 

2.006 

1169.16 

12.06 

4213.94 

114.02 

50.86 

9.56 

63.40 

8.66816 

2.539 

1175.18 

11.30 

4759.95 

146.24 

59.18 

12.15 

63.74 

8.67094 

2.600 

1175.38 

-1.67 

4823 .12 

149.88 

59.18 

12.45 

63.75 

8.67126 

3.112 

1173.38 

4.36 

5349 . 12 

180 . 16 

59.09 

14.94 

63 . 66 

8.67393 

3.645 

1164.48 

-19.20 

5893.98 

211.51 

58.65 

17.52 

63.22 

8.67671 

4.054 

1153.28 

-33.32 

6309 . 70 

235.44 

58.08 

19 . 50 

62 . 66 

8 . 67885 

4 . 382 

1142.12 

-33.80 

6639.00 

254 .38 

57.54 

21.00 

62.09 

8 . 60056 

4.648 

1130.45 

-48.58 

6904.06 

269 . 62 

56.94 

22.35 

61.50 

8 . 68195 

4.894 

1118.50 

-40.41 

7146.19 

283.54 

56 . 35 

23.51 

60.08 

8.68324 

5.099 

1108.47 

-  55 . 57 

7346.10 

295.03 

55.84 

24.47 

60.37 

0 . 68430 

5.283 

1097.31 

-61.91 

7524 . 33 

305.28 

55.28 

25.33 

59.80 

8.68527 

5.447 

1087.16 

-61.01 

7681.28 

314.29 

54.78 

26.10 

59.27 

8.60612 

5.611 

1077.07 

-61.42 

7836.83 

323.23 

54.28 

26.05 

50.75 

8.68698 

5.774 

1066.62 

-70.73 

7990.99 

332.08 

53.75 

27.61 

58.21 

8.68784 

5.930 

1054.50 

-74.37 

8143 . 54 

340.03 

53.14 

28.35 

57.59 

8.60069 

6.102 

1042.41 

-74.03 

0294.40 

349.49 

52.54 

29.09 

56.96 

8.68955 
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Table  3s  (cont.) 
ROCKET  SAMPLE  OUTPUT 


TIME 

PC 

FOOT 

FREE  VOL 

6.266 

1030.32 

-73.41 

8443.60 

6.360 

1019.22 

-25894.06 

8527.95 

6.360 

1006.86 

-25580.20 

0527.95 

6.361 

990.63 

-25167.67 

8527.95 

6.361 

971,65 

-21761.79 

8527.95 

6.362 

958.93 

-24362.40 

8527.95 

6.363 

928.25 

-23582.91 

8527.95 

6.365 

898.54 

-22828.35 

8527.95 

6.366 

869.79 

-22097.94 

8527.95 

6.368 

815.03 

-20706.56 

0527.95 

6.371 

763.71 

-19102.71 

0527.95 

6.376 

670.56 

-17036.25 

8527.95 

6.381 

588.77 

-14958.36 

8527.95 

6.386 

516.95 

-13133.84 

8527.95 

6.391 

153.90 

-11531.05 

8527.95 

6.397 

398.53 

-10125.25 

8527.95 

6.102 

319.92 

■0890.23 

0527.95 

6.107 

307.21 

-7805.81 

8527.95 

6.112 

269.76 

-6853.72 

8527.95 

6.117 

236.86 

-6017.73 

8527.95 

6.122 

207.96 

-5283.71 

8527.95 

6.127 

182.60 

-4639.22 

8527.95 

6.132 

160.32 

-1073.35 

8527.95 

6.138 

140.77 

-3576.49 

8527.95 

6.113 

123.60 

-3110.21 

8527.95 

6.118 

100.52 

-2757.20 

8527.95 

6.153 

95.20 

-2420.88 

8527.95 

6.458 

83.66 

-2125.58 

8527.95 

6.163 

73.16 

-1866.31 

8527.95 

6.173 

56.63 

-1438.77 

0527.95 

6.184 

43.66 

-1109.18 

8527.95 

6.191 

33.65 

-855.08 

8527.95 

6.509 

22.78 

-578.79 

8527,95 

AVERAGE 

PRESSURE 

FROM 

000  SEC. 

MAXIMUM 

PRESSURE 

WAS  1175.41  F.SI 

TIME 

PC 

PE 

F 

.000 

11 . 70 

.  00 

.  00 

.001 

27.75 

.00 

.  00 

.001 

30.26 

5.18 

273.95 

.001 

51.96 

1 . 87 

417.95 

.002 

66.64 

1.63 

577.29 

.002 

76.89 

1.50 

690.08 

.002 

87.46 

4.39 

809.61 

.  003 

98.30 

1 . 29 

933.04 

.003 

109.38 

4.20 

1060.59 

.003 

120.68 

4.11 

1191.83 

.001 

132.16 

4.01 

1326.36 

.001 

113.80 

3.97 

1163.80 

P  WBIGtIT 

WDOTP 

INS  WT 

WDOT 

AT 

350 

.05 

51.94 

29. 

83 

56, 

31 

8 

69040 

362 

.89 

.  00 

30. 

25 

55. 

75 

8 

69009 

362 

.89 

.00 

30. 

25 

55. 

07 

8. 

69090 

362 

.89 

.00 

30. 

25 

51. 

19 

8 

69090 

362 

.89 

.00 

30. 

25 

S3 

31 

8 

69090 

362 

.89 

.00 

30. 

25 

52. 

IS 

0. 

69091 

362 

.89 

.00 

30. 

25 

50. 

77 

8. 

69092 

362 

.89 

.  00 

30. 

25 

19  . 

15 

8. 

69092 

362 

.89 

.00 

30. 

25 

17. 

58 

8 

69093 

362 

.89 

.00 

30. 

25 

11. 

58 

0. 

69091 

362 

.89 

.00 

30. 

25 

41. 

77 

8 

69096 

362 

.89 

.00 

30. 

25 

36. 

68 

8. 

69090 

362 

.89 

.00 

30. 

25 

32. 

21 

8. 

69101 

362 

.09 

.00 

30. 

25 

28. 

28 

8. 

69101 

362 

.89 

.00 

30. 

25 

21  . 

03 

8  . 

69106 

362 

.89 

.00 

30. 

25 

21  . 

80 

8. 

69109 

362 

.89 

.00 

30. 

25 

19  . 

11 

8. 

69112 

362 

.89 

.00 

30. 

25 

16. 

81 

8. 

691 11 

362 

.89 

.00 

30. 

25 

11  . 

76 

8. 

69117 

362 

.09 

.00 

30. 

25 

12. 

96 

8. 

69120 

362 

.  89 

.00 

30. 

25 

11. 

38 

8. 

69122 

362 

.89 

.00 

30. 

25 

9  . 

99 

8 

69125 

362 

.09 

.00 

30. 

25 

8. 

77 

8. 

69128 

362 

.89 

.00 

30. 

25 

7. 

70 

8. 

69130 

362 

.89 

.00 

30. 

25 

6. 

7(5 

8. 

69133 

362 

,89 

.00 

30. 

25 

5. 

91 

8, 

69136 

362 

.89 

.00 

30. 

25 

5. 

21 

0. 

69138 

362 

.89 

.00 

30. 

25 

1 . 

58 

0. 

69111 

362 

.89 

.00 

30. 

25 

1 . 

02 

8. 

69141 

362 

.89 

.  00 

30. 

25 

3  . 

10 

8. 

69119 

362 

.89 

.00 

30. 

25 

2. 

39 

8. 

69151 

362 

89 

.00 

30. 

25 

1  . 

81 

8. 

69160 

362 

89 

.00 

30. 

25 

1 . 

25 

0. 

69168 

TO 

6.509 

SEC. 

WAS 

1110.9 

M3 

r.si 

WRIGHT 

CF 

WUOT  EX 

ISP 

.00 

.000 

.00 

.00 

.01 

.000 

.00 

.  00 

.02 

.827 

.00 

113.31 

.03 

.929 

.00 

120.11 

.01 

1.001 

.00 

127.75 

.05 

1.038 

.00 

132.32 

.06 

1 . 069 

.01 

136.51 

.07 

1.096 

.01 

110.11 

.08 

1 .120 

.01 

113.99 

.09 

1.111 

.01 

117.30 

.  10 

1.159 

.02 

150.38 

.11 

1.176 

.02 

153.21 
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Table  3;  (cont.) 
ROCKET  SAMPLE  OUTPUT 


TIME 

PC 

PE 

F 

.004 

155.57 

3.91 

1603.01 

.005 

167.45 

3.85 

1746.08 

.005 

191.48 

3.75 

2036.21 

.006 

215.74 

3.66 

2353.98 

.007 

210 . 11 

3.59 

2676 . 77 

.007 

261.19 

3.52 

3002 . 53 

.000 

200 . 00 

3.16 

3329 . 71 

.009 

312.95 

3.56 

3674.16 

.009 

336.09 

3.83 

1033 . 64 

.010 

360.55 

4 . 10 

4389.03 

.010 

383.89 

4.37 

1739.64 

.011 

406.88 

4.63 

5084 . 88 

.012 

451.65 

5.14 

5757.28 

.014 

494.66 

5.63 

6403.28 

.015 

535.78 

6.10 

7020.91 

.016 

571 .91 

6.54 

7609.00 

.010 

612.09 

6.97 

8167.05 

.  019 

617.24 

7.37 

8695.02 

.020 

680.42 

7 . 74 

9193.27 

.021 

711.65 

0.10 

9662.43 

.023 

741.01 

0.43 

10103.34 

.024 

768.55 

0.75 

10517.02 

.025 

794 . 35 

9.04 

10904.56 

.026 

018.49 

9.31 

11267.14 

.029 

862.10 

9.81 

11922.16 

.032 

900.05 

10.24 

12492.06 

.034 

932.96 

10.62 

12986.34 

.037 

961.43 

10.94 

13414.09 

.039 

986.03 

11.22 

13783.48 

.  042 

1007.21 

11.46 

11102.08 

.014 

1025.51 

11.67 

11376 . 50 

.047 

1041.23 

11.85 

11612.75 

.050 

1054.76 

12.00 

14815.98 

.052 

1066.40 

12.14 

11990.82 

.055 

1076.41 

12.25 

15141 . 17 

.060 

1092.13 

12.43 

15381.87 

.065 

1104.32 

12.57 

15560.39 

.073 

1116.71 

12.71 

15746.61 

.083 

1126.81 

12.82 

15898.38 

.  106 

1137.17 

12.94 

16054.16 

.183 

1148 . 30 

13.07 

16222 . 20 

.  249 

1155.79 

13 . 16 

16335.24 

.255 

1110.96 

12.99 

16112.51 

.263 

1128.01 

12.84 

15918.07 

.273 

1117.77 

12 . 72 

15764.10 

.  311 

1107.56 

12.61 

15611.14 

.  322 

1107.11 

12.61 

15609.22 

.532 

1118.03 

12.71 

15770.61 

WEIGirr 

CF 

WDOT  EX 

ISP 

.12 

1 

.  191 

.02 

155 

.91 

.13 

1 

.201 

.02 

158 

.41 

.15 

1 

.228 

.03 

162 

.96 

.  18 

1 

.  260 

.01 

167 

.  19 

.20 

1 

.  200 

.05 

1  71 

.  21 

.23 

1 

.  311 

.06 

175 

.  05 

.25 

1 

.  332 

.  07 

170 

.  61 

.  28 

1 

.356 

.08 

181 

.99 

.31 

1 

.  303 

.09 

185 

.31 

.33 

1 

.406 

.10 

188 

.54 

.36 

1 

.126 

.12 

191 

.66 

.  39 

1 

.  143 

.13 

194 

.63 

.  45 

1 

.472 

.16 

200 

.  10 

.51 

1 

.195 

.20 

204 

.95 

.57 

1 

.514 

.23 

209 

.  21 

.63 

1 

.529 

.  27 

213 

.  03 

.  70 

1 

.  541 

.31 

216 

.10 

.  76 

1 

.552 

.  36 

219 

.  10 

.  83 

1 

.561 

.11 

222 

.08 

.90 

1 

.  568 

.46 

224 

.  19 

.97 

1 

.575 

.51 

226 

.66 

1 

.04 

1 

.501 

.56 

228 

.  62 

1 

.11 

1 

.506 

.62 

230 

.  10 

1 

.  18 

1 

.590 

.67 

232 

.02 

1 

.  32 

1 

.597 

.  79 

231, 

.  86 

1 

.47 

1 

.603 

.91 

237, 

.25 

1 

.62 

1 

.608 

1 

.04 

239 

.29 

1 

.76 

1 

.611 

1 

.18 

241 

05 

1 

.92 

1 

.615 

1 

.31 

212  . 

57 

2 

.07 

1 . 

617 

1 

.45 

213  . 

90 

2 

.22 

1 

619 

1 

.  59 

215. 

07 

2 

.37 

1 . 

621 

1 

.71 

216. 

10 

2 

.53 

1 . 

622 

1 

.88 

217. 

02 

2 

.68 

1  . 

624 

2 

.03 

217. 

84 

2 

.84 

1 . 

625 

2 

.  18 

218  . 

57 

3 

.  15 

1  . 

626 

2 

.48 

249. 

83 

3 

.46 

1 . 

627 

2 

.  79 

250. 

87 

3 

.94 

1  . 

629 

3 

.25 

252  . 

12 

4 

.57 

1  . 

630 

3 

.07 

253  . 

37 

6 

on 

1 

C  1  1 

5 

.29 

255. 

19 

10. 

79 

1  . 

632 

10 

.05 

257. 

68 

11. 

97 

1  . 

632 

11  . 

21 

258. 

55 

15. 

33 

1  . 

631 

14  . 

57 

258. 

60 

15. 

80 

1  . 

630 

15 

05 

258  . 

66 

16. 

12 

1  . 

629 

15. 

67 

258. 

71 

18. 

71 

1  . 

628 

17. 

99 

258. 

83 

19. 

35 

1 . 

628 

18. 

60 

258. 

86 

32. 

10 

1. 

629 

11  . 

29 

259. 

19 

30 


Table  3!  (cont.) 
ROCKET  SAMPLE  OUTPUT 


TIME 

PC 

.736 

1128.33 

.982 

1138.11 

1.228 

1148.53 

1.576 

1158.81 

2.006 

1169 .16 

2.539 

1175.18 

2.600 

1175.38 

3.112 

1173.38 

3.615 

1164.18 

4.051 

1153.28 

1.382 

1142.12 

1 . 618 

1130.15 

1.891 

1118.50 

5.099 

1 108 . 17 

5.283 

1097.31 

5. 147 

1007.16 

5.611 

1077.07 

5.771 

1066 . 62 

5.938 

1051.58 

6.102 

1012.11 

6.266 

1030.32 

6.360 

1019.22 

6.360 

1006 . 86 

6.361 

990 . 63 

6.361 

971.65 

6.362 

958.93 

6.363 

928.25 

6.365 

898.51 

6.366 

069 . 79 

6.368 

815.03 

6.371 

763.71 

6.376 

670.56 

6.381 

588 . 77 

6.306 

516.95 

6.391 

153.90 

6.397 

398.53 

6.102 

319 . 92 

6 . 107 

307.21 

6.112 

269 . 76 

6.117 

236.86 

6 . 122 

207.96 

6 . 127 

182  60 

6.132 

160.32 

6.138 

110.77 

6 . 113 

123.60 

6.118 

108.52 

6.153 

95.28 

6.150 

83 .66 

PE 

F 

12.85 

15927.39 

12.97 

16081.13 

13.08 

16235.53 

13.20 

16393.98 

13.33 

16553.22 

13.10 

16619.05 

13.40 

16652.59 

13.39 

16627.51 

13.29 

16198.77 

13.17 

16331.23 

13.01 

16169.36 

12.91 

15996.25 

12.78 

15818.51 

12.66 

15669.13 

12.51 

15503.05 

12.12 

15351 . 63 

12.31 

15200.99 

12.19 

15015.10 

12.06 

11865.09 

11.92 

11683 .13 

11.78 

11502.13 

11.66 

11335.81 

11  .51 

11119.61 

11.33 

13901 .90 

11.15 

13661 . 10 

10 .9  7 

13127.16 

10.62 

12961 . 72 

10.28 

12517.07 

9.95 

12083.71 

9.32 

11250.29 

8.73 

10181.76 

7.67 

9080.83 

6.73 

7818.10 

5.91 

6765.68 

5.19 

5815.28 

1.56 

1980.80 

1.00 

1210 . 10 

3.51 

3601 . 78 

3.50 

3085.91 

3.60 

2611 .11 

3.69 

2261 .16 

3 . 79 

1936.15 

3.89 

1667.17 

3.99 

1133.56 

1 . 09 

1230.77 

1 . 20 

1051 . 80 

1 . 31 

902.10 

1.13 

769 . 88 

WEIGHT 

CF 

41 

.65 

1 

630 

59 

.81 

1 

631 

75 

.17 

1 

632 

97 

.08 

1 

632 

121 

.37 

1 

633 

158 

.39 

1 

631 

162 

.32 

1 

631 

195 

.09 

1 

631 

229 

.01 

1 

633 

251 

.91 

1 

632 

275 

.15 

1 

631 

291 

.97 

1 

630 

307 

.05 

1 

629 

319 

.51 

1 

628 

330. 

.61 

1 

627 

310 

39 

1 

626 

350  . 

08 

1 

625 

359  . 

68 

1  . 

621 

369. 

19 

3  . 

622 

378. 

59 

1  . 

621 

38  7. 

98 

3  . 

620 

393. 

11 

1  . 

638 

393. 

1 1 

3  . 

637 

393. 

11 

1  . 

615 

393. 

1 1 

3  . 

613 

393. 

11 

3  . 

611 

393. 

11 

1  . 

607 

393. 

11 

1  . 

603 

393. 

11 

3  . 

599 

393. 

11 

1  . 

589 

393. 

11 

1  . 

580 

393. 

11 

1  . 

558 

393 

.  11 

1 

.531 

393 

11 

1 

.506 

393 

11 

1 

.  171 

393 

11 

1 

.138 

393 

11 

1 

.39  7 

393 

11 

1 

.350 

393 

11 

1 

.316 

393. 

11 

1 

.205 

393. 

11 

1 

.251 

393. 

11 

1  . 

220 

393  . 

11 

1  . 

.196 

393. 

11 

1  . 

172 

393. 

11 

1  . 

116 

393. 

11 

1  . 

118 

393. 

11 

1  . 

009 

393. 

11 

1  . 

059 

WDOT  EX 

ISP 

13 

.  77 

259 

.11 

58 

.  09 

259 

.61 

74 

.11 

259 

.  79 

95 

.92 

259 

.99 

123 

.08 

260 

.21 

156 

.92 

260 

.11 

160 

.81 

260 

.13 

193 

.46 

260 

.56 

227 

.21 

260 

.63 

253 

.02 

260 

.65 

273 

.16 

260 

.  65 

289 

.  91 

260 

.  62 

301 

.  95 

260 

.  59 

317 

.  37 

260 

.  56 

328. 

11 

260. 

.52 

338. 

20 

260. 

18 

317. 

86 

260  . 

13 

357. 

15 

260. 

38 

366. 

93 

260. 

33 

3  76. 

32 

260. 

27 

385  . 

60 

260. 

21 

390. 

05 

260. 

17 

390. 

87 

260. 

17 

390. 

91 

260. 

1  7 

390. 

91 

260. 

17 

390. 

9  7 

260. 

17 

391. 

01 

260. 

16 

391. 

10 

260. 

15 

391  . 

17 

260. 

16 

391 . 

28 

260. 

16 

391  . 

10 

260. 

16 

391.60  260. J5 
391.77  260. IS 
391.93  260.11 
392.06  260.13 
392.10  260.12 
392.29  260.11 
392.38  260.10 
392.16  260.09 
392.53  260.00 
392.59  260.07 
392.65  260.07 

392.70  260.06 

392.71  260.05 
392.70  260.01 
392.81  260.01 
392.81  260.03 
392.06  260.02 
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Table  3:  (cent.) 
ROCKET  SAMPLE  OUTPUT 


TIME 

PC 

PE 

F 

WEIGHT 

CF 

WDOT  EX 

ISP 

6.463 

73.46 

4.54 

655.25 

393.14  1 

.  026 

392.88 

260.02 

6.473 

56.63 

4 . 79 

470.06 

393 . 14 

.955 

392.92 

260.01 

6.484 

43.66 

5 . 04 

331.44 

393.14 

.874 

392.95 

260.00 

6.494 

33 . 65 

5.31 

227.84 

393.14 

.779 

392.97 

259.99 

6.509 

22.78 

5.74 

119.40 

393 . 14 

.603 

392.99 

259.99 

AVERAGE 

THRUST  FROM 

000 

SEC.  TO 

6 . 509 

SEC. 

WAS  15696.570  Ibf 

TOTAL  IMPULSE= 

102172 

.  50 

SEC. 

THRUST  VS 

!.  TIME 

THRUST 

(1000 

Ibf ) 

2  . 

4  . 

6 

8  . 

10.  12. 

14 

16  . 

18.  20. 

.000 
.052 
.106 
.  183 
.249 
.311 
.  532 
.  736 
.982 
1.228 
1.576 
2.006 
2.539 
2.600 
3 .112 
3.645 
4.054 
4 . 382 
4 . 648 
4 . 894 
5.099 
5.283 
5.447 
5.611 
5.774 
'=.938 

6.  x02 
6.266 
6.360 
6.412 
6.463 
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Table  3s  (cont.) 
ROCKET  SAMPLE  OUTPUT 


rRKSSURF  VS.  TTMR 

TRESSURE 
(1000  ppi) 


.000 
.052 
.106 
.193 
.249 
.311 
.532 
.736 
.982 
1.228 
1.576 
2.006 
2.539 
2.600 
3.112 
3.645 
4.054 
4  ,  392 
4.649 
4.894 
5.099 
5.283 
5.44  7 
5.611 
5.7  74 
5.938 
6.102 
6.266 
6.360 
6.412 
6.463 
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IV.  RESULTS  AND  CONCLUSIONS 


A.  RESULTS 

In  order  to  verify  the  accuracy  of  the  modifications 
made  to  NEWPEP,  several  problems  were  run  using  the  old 
version,  the  updated  version  and  the  NASA  SP-273  code  [Ref- 
6] .  The  results  for  the  following  two  problems  are 
summarized  in  Table  4. 

1.  A  mixture  of  6%  by  weight  H2  and  94%  by  weight  O2 
was  burned  at  a  chcimber  pressure  of  3.447  MPa  (500  psi)  and 
expanded  to  an  eimbient  pressure  of  0.1014  MPa  (14.7  psi). 

2 .  A  mixture  of  85%  by  weight  ammonium  perchlorate  and 
15%  HTPB  (Sinclair)  was  burned  at  6.8947  MPa  (1000  psi)  and 
expanded  to  an  ambient  pressure  of  0.1014  MPa  (14.7  psi) . 
There  were  no  significant  differences  in  any  of  the  output 
values . 

B.  CONCLUSIONS 

PEP93  is  an  interactive  program  that  works  well  on  a 
personal  computer.  It  allows  the  user  to  evaluate  a  wider 
variety  of  performance/design  problems  while  using  an 
existing,  familiar  program. 

ROCKET  was  successfully  modified  to  allow  the  use  of 
multiple  grains.  This  greatly  enhances  the  ability  to 
predict  the  performance  of  actual  rocket  motors. 
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Table  4: 

SAMPLE  COMPARISON  OF  NEWPEP/PEP93/NASA  SP-273 


Sample  I 


Ingredients 


H2  6%  by  weight 


02  94%  by  weight 


NEWPEP 


Chamber 


P  (MPa) 


T(K| 


HIkJ/kgl 


SIkJ/kgKI 


Exhaust  P  {MPal 


TIKI 


HIkJ/kg) 


SIkJ/kgKI 


Performance 


Frozen 


Specific 


Impulse  (sec) 


Im/sec) 


Ae/A* 


Cf 


Shifting 


Specific 


Impulse  IsecI 


Im/secl 


Ae/A* 


Cf 


ShiftingIFrozen 


Spedfic 


Impulse  IsecI 


|m/sec| 


3.447 

3.447 

3271.9 

3273.05 

12.268 


1.989 


30961 


■675.298 


12.268 


2192 


3636.31 


12.268 


2.268 

'i.99 


3096.5 


675.501 


12.268 


2191.9 


-3636.241 


12.368 


12.2631 


1.989 


12J631 
0. 


2191.79 


3637.1 


12.2631 


275 


1772.7 


5.77 


274.9 


1773 


5.7662 

1.521 


Sample  2 


Ingredients 


AP  85%  by  weight 


HTPB  15%  by  weight 


NEWPEP  PEP93  NASA 


SP273 


6.895 


2844 


2132.79 


10.15 


3.885 


2595.8 


2133.15 


10.146 


3.885 


2596.08 


1324 

1323.8 

1323.67 

-4995.7 

'4995.56 

-4996.42 

10.15 

10.15 

10.146 

241.1 

241 

~  '1501 

1501 

8.091 

8.0894 

1.575 

1.575 

243.9 


1514 


8.3596 


1.581 


APPENDIX  A 


USER'S  GUIDE  FOR  PEP93 

(Changes  in  Bold) 

A.  GENERAL 

1.  The  files  NEWPEP.FOR,  NEWEQUIL.FOR  and  NEVTOSIGN.FOR 
are  needed  to  create  NEWPEP.EXE.  In  addition,  several 
common  blocks  are  needed.  BLANK. INC,  CHARA.INC,  IBRIUM.INC, 
SCRATCH. INC  and  SIZE. INC. 

2.  The  data  files  JANNAF.DAT  and  PEPCODED.DAT  are 
needed  to  run  PEP93. 

B.  INPUT 

1.  Up  to  15  ingredients  may  be  input  by  the  user. 

These  ingredients  may  be  coded  ingredients  from  the 
PEPCODED.DAT  file  or  user  defined  ingredients,  or  both. 

2.  The  user  is  prompted  for  the  following  inputs  when 
PEP93  is  run: 

a.  File  naime  (up  to  30  characters) 

b.  Run  title  (up  to  17  characters) 

c.  Vitiated  air  option 

(1)  Number  of  coded  Ingredients  for  vltlator 

(2)  Number  of  user  defined  ingredients  for 
vltlator 

(3)  Ingredient  code  numbers  from  PEPCODED.DAT 
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(4)  User  defined  Ingredient  name,  con^osition, 
heat  of  formation  and  density 

(5)  Chamber  pressure 

(6)  Chamber  inlet  temperature 

(7)  Ingredient  weights 

d.  Number  of  coded  ingredients  for  combustor 

e.  Number  of  user  defined  ingredients  for  combustor 

f.  Niimber  of  runs 

g.  Density  exponent  used  in  calculating  density 
(defaults  to  1.0) 

h.  Propellant  temperature  to  adjust  the  system 
enthalpy  for  heating  or  cooling 

(default  is  298K) 

i.  Quadratic  coefficients  for  solid  specific  heat. 
Used  to  adjust  the  system  enthalpy  for  heating 
or  cooling  (defaults  to  0.3  cal/gm-K) 

j.  Nozzle  discharge  coefficient  (defaults  to  1.0) 

k.  Options 

(1)  Delete  exit  calculations 

(2)  Include  ionic  species  in  calculations 

(3)  Include  boost  velocities  and  nozzle  design 
data 

(4)  Input  pressures  in  psi  Instead  of  MPa 

(5)  Increase  numerical  precision  of  species 
concentrations 
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(6)  Output  a  list  of  all  combustion  species 
considered 

(7)  Fix  temperature 

(8)  Debug  options 

(a)  Thermo  data  at  every  guess 

(b)  Values  of  J,  M,  VF,  PR,  VA  in 
subroutine  TWITCH 

(c)  Species  composition  every  guess 

(d)  Log  of  equilibrium  constant  every 
guess 

(e)  Classification  of  species  each 
iteration  from  TWITCH 

(9)  Exact  or  approximate  throat  calculations 

(10)  Normalize  or  not  normalize  to  100  grains 

(11)  Calculate  non- ideal  expansion  performance 

(12)  Calculate  chamber  Hach  number 

l.  Ingredient  code  numbers  from  PEPCODED.DAT 

m.  User  defined  ingredient  name,  composition,  heat 
of  formation  and  density 

n.  Chamber  pressure (MPa) 

o.  Exhaust  pressure  or  area  ratio  and  ambient 
pressure  (see  note  2) 

p.  Ingredient  weights 

q.  Ratio  of  chamber  area  to  throat  area 
3.  Notes: 

a.  Format  for  user  defined  ingredients 
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C.  SAMPLE  INPUT 

TYPE  IN  OUTPUT  FILE  NAME  (DEFAULT:  PEPOUT.DAT) 

If  you  wish  to  direct  to  a  printer,  type  in 
printer  name  such  as  "LPTl"  for  the  IBM  PC 
or  "LPAO:"  for  the  VAX. 

To  send  output  to  screen,  type  "TT; "  (VAX  &  IBM  PC) 


TYPE  RUN  TITLE  (char* 17)  OR  type  quit, 
sample 

DO  YOU  WANT  TO  CALCULATE  COMPOSITION  AND  ENTHALPY  OF 
VITIATED  AIR?  l*YES,2=NO 
1 

INPUT  NO.  OF  CODED  ING.,  NO.  OF  USER  DEFINED  ING.  FOR 
VITIATOR 
2,0 

READ  IN  2  INGREDIENT  CODE  NUMBERS  SEP.  BY  COMMAS 
44,729 

READ  IN  CHAMBER  PRESSURE  (MPa) ,  CHAMBER  INLET  TEMPERATURE 
(K) ,  wtl,wt2,  +etc.  2.068,700. ,96. ,4. 

BEGIN  ENGINE  DATA  INPUT. 

NO.  OF  CODED  ING.,  NO.  OF  USER  DEFINED  ING.,  AND  NO. OF  RUNS 

ss 

NOTE:  INPUT  VITIATED  AIR  AS  A  USER  DEFINED  INGREDIENT. 

1,1,1 

DENSITY  EXPONENT  (DEFAULTS  TO  1.0)  = 

PROPELLANT  TEMPERATURE  (DEFAULTS  TO  298.0)  = 

ENTER  THREE  VALUES  FOR  CSUBP  FIT 
FIRST  VALUE  =0.0  DEFAULTS  TO  0.30 

NOZZLE  DISCHARGE  COEFFICIENT  (DEFAX7LTS  TO  1.0)  = 

OPTIONS;  1- DELETE  EXIT  CALCULATIONS 

2 - INCLUDE  IONIC  SPECIES  IN  CALCULATIONS 

3 - INCLUDE  BOOST  VELOCITIES  AND  NOZZLE  DESIGN  DATA 

4 - INPUT  PRESSURES  IN  PSI  INSTEAD  OF  MPA 

5- "N”  ORDERS  OF  MAGNITUDE  SPECIES  PRECISION 
INCREASE 

6 - OUTPUT  A  LIST  OF  ALL  COMBUSTION  SPECIES 
CONSIDERED 

7 - DEBUG  OPTIONS 

8- "P-T-H-S"  MAP  OPTION 

9- "0"=exact  throat  "l"=approx  throat  "2"=exact 
with  throat  composit . 

0-"l"=do  not  normalize  weights  to  100  gms 
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"O"=nonnalize 

A-CALCULATE  NON- IDEAL  EXPANSION  PERFORMANCE 
B- CALCULATE  CHAMBER  MACH  NUMBER 

1234567890AB**put  0-no  or  1-yes  under  number  (opt  5  takes  0 
or"n") 

000000000011 

READ  IN  1  INGREDIENT  CODE  NUMBERS  SEP.  BY  COMMAS 
714 

VITIATED  AIR  IS  AUTOMATICALLY  ENTERED  AS  FIRST  USER  DEFINED 
INGREDIENT 

READ  IN  CHAMBER  PRESSURE  (MPa) , AREA  RATIO, AMBIENT  PRESSURE 
(MPa),  WT1,WT2,  +ETC. 

INCLUDE  DECIMAL  POINT  AND  SEPARATE  BY  COMMAS 

INPUT  A  ZERO  CHAMBER  PRESSURE  TO  ADD  INGREDIENTS  AGAIN 

INPUT  A  -1.0  TO  END  PROGRAM 

1.379,5. ,14.7,7. ,93. 

INPUT  RATIO  OF  CHAMBER  AREA  TO  THROAT  AREA  TO. 

CALCULATE  MACH  NUMBER  AT  STATION  4 
1.5 

TYPE  RUN  TITLE  (char* 17)  OR  type  quit 
quit 

Stop  -  Program  terminated 


D .  SAMPLE  OUTPUT 

See  Table  2 . 


APPENDIX  B 


USER'S  GUIDE  FOR  ROCKET 

A.  INPtJT 

1.  The  data  file  CARD.DAT  is  needed  to  run  ROCKET  and 


consists 

of 

the  following  variables. 

Title 

Ath 

Ae 

Vol  Tp 

Po 

Pb 

Pi 

Vcf 

y  Pd 

td 

ET# 

pressure  1 

erosion  rate 

at 

pressure 

1 

pressure  2 

erosion  rate 

at 

pressure 

2 

etc- 

etc. 

G# 

- 

« 

r 

n 

Pref 

c* 

X 

Pp 

ti 

^ab 

pts 

id# 

Wb 

Ab 

0  0 

0 

where 

Ath  Throat 

area  (in^) 

Enter 

a  positive  area 

if 

it  is  known 

Enter  a  negative  area  for  a  first  guess  if  unknown 
Enter  zero  if  you  have  no  first  guess  (very  slow) 
Ae  Nozzle  exit  area  (in^) 
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Vol  Total  motor  voliime  without  propellant  or  igniter  (in^) 
Tp  Ambient  propellant  temperature  (°F) 

Po  Ambient  pressure  (psia) 

Pb  Throat  plug  closure  blowout  pressure  (psia) {typically 
35} 

Pi  Initial  motor  pressure  (psia) 

Nozzle  loss  coefficient 

7  Generally  use  7^  for  shifting  equilibriiam  flow 
(PEP93) 

Pd  Design  pressure  for  sizing  the  throat  (psia) 
td  Throat  erosion  delay  time  (sec) 

Positive  value;  fixed  rate  from  t=0  (i.e.  0.0005 
in/sec) 

Negative  :  equals  the  time  delay  for  onset  of  erosion 
Negative  or  0:  must  include  pressure  vs.  erosion  rate 
table  (Erosion  rates  are  radial  in  in. /sec) 

ET#  Nximber  of  points  for  throat  erosion  table 
G#  Number  of  grains 

r  Propellant  burning  rate  at  P^-^f  (in. /sec) 
n  Burning  rate  exponent  (r=aPj,“) 

IIj^,  Propellant  temperature  sensitivity  (%/°F) 

Pj.gf  Reference  pressure  for  the  specified  burning  rate 
(psia) 

C*  Characteristic  exhaust  velocity  at  P^ef  (ft/sec) 

X  Pressure  correction  exponent  for  C* 
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Pp  Density  of  grain (Ibm/in^) {propellant, igniter  or 
insulation} 

t^  Ignition  time  (sec)  {O  if  no  igniter} 

At^  Ignition  time  delay  (sec)  {0  if  no  igniter} 
pts  Number  of  points  in  web  burned  vs.  burning  surface 
area  table 

Wb  Web  distance  burned  (in) 

Ab  Propellant  burning  surface  area  (in^) 

Ab^  Exposed  insulation  area  (in^) 
id#  Identification  number  for  grain 
1 : igniter 
2 : insulation 
3 : propellant 

rgj3  Insulation  ablation  rate 
2 .  Notes : 

a .  All  input  data  is  formatted  F12 . 6  except : 

(1)  Title  -  no  format  required 

(2)  ET#,G#  -  format  13 

b.  Enter  0  in  columns  3,6  and  9  of  the  last  row. 

c.  Igniters: 

(1)  r  always  equal  to  1 

(2)  pts  is  #  of  points  in  time  vs.  m  curve 

(3)  Wb=time  (sec) 

(4)  Ab=  m  (Ibm/sec) 

d.  Insulation: 

(1)  r,n,n]^  and  Pref  refer  to  the  propellant 
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grain,  not  to  the  insulation.  They  are  used  to  calculate  a 
reference  web  which  can  be  used  to  describe  the  amount  of 
insulation  exposed  during  a  firing.  C*,X,Pp  ,  t and  At^^  are 
for  the  insulation. 

(2)  Wb  is  the  web  of  the  associated  propellant 
grain.  Ab  is  the  exposed  insulation  area. 

B .  OUTPUT 

1.  Input  data  is  printed  with  all  values  labeled. 

2.  Output  data  consists  of  the  following: 

a.  Total  propellant  weight  and  igniter  weight. 

b.  A  table  of  time,  chamber  pressure,  change  in 
pressure,  total  motor  volume,  propellant  weight  consumed, 
rate  of  mass  addition  due  to  propellant  burning,  insulation 
weight  cons\imed,  rate  of  mass  flow  through  the  nozzle  and 
throat  area. 

c.  Average  and  maximxim  pressure  for  run. 

d.  A  table  of  time,  chamber  pressure,  exit 
pressure,  thrust,  coefficient  of  thrust,  rate  of  mass  flow 
through  the  nozzle  and  specific  impulse. 

e.  Average  thrust  and  total  Inpulse. 

f.  Plots  of  thrust  versus  time  and  pressure  versus 
time. 

3.  Time,  pressure,  and  thrust  are  also  put  into  the 
data  file  PLOT.DAT 
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C.  SAMPLE  INPUT 


See  Table  5 . 

D.  SAMPLE  OUTPUT 

See  Table  3 . 
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Table  5s 

ROCKET  SAMPLE  INPUT 


iSpercent 

e. 657700  105.450000  8042.477000  70.000000  10.130000  35.000000 


14.700000  0.960600 

3 

0.580000  0.350000 

0. 062300  0.000000 

11  1 

0.000000  1448.300000 
0.380900  1519.500000 
0.761800  1542.500000 
1.142800  1557.300000 
1.523700  1563.900000 
1.904600  1562.300000 
2.285500  1552.500000 
2.666400  1534.500000 
3.047400  1508.200000 
3.428300  1473.800000 
3.809200  1431.200000 
0.580000  0.350000 

0.062300  0.000000 

tl  2 

0.000000  1488.300000 
0.380900  1519.500000 
0.761800  1542.500000 
1.142800  1557.300000 
1.523700  1563.900000 
1.904600  1562.300000 
2.285500  1552.500000 
2.666400  1534.500000 
3.047400  1508.200000 
3.428300  1473.800000 
3.809200  1431.200000 
1.000000  0.000000 
0.062300  0.000000 

2  1 

0.000000  1.400000 

0.250000  1.900000 

0  0  0 


1.144300  10.30.000000 

0.200000  1030.000000 
0.000000 


0.200000  1030.000000 
0.000000  0.050000 


o.ooooon  0.000000 


0 . 000500 

5106.660000  0.045000 


5106.660000  0.045000 


5106.660000  0.045000 
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APPENDIX  C 


PEP93  SUBROirriNE  EQUIVALENCE 


SUBROUTINE  EQUIVALENCE (equiv) 
INCLUDE  'size.inc' 

INCLUDE  'blank. inc' 

INCLUDE  ' chara . inc ' 

INCLUDE  ' scratch . inc ' 
INCLUDE  ' ibrium . inc ' 
real  ostoic, oin, equiv 
coiranon/llist/atwt  (100) 
ostoic=0 . 0 


do  10  j=l,is 

if (aspec ( j ) 
if (aspec  ( j ) 
if (aspec (j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec  { j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec  ( j ) 
if  (aspec  (j ) 
if (aspec ( j ) 
if (aspec (j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec  ( j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if  (aspec  (j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec  (j ) 
if (aspec  ( j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec  ( j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec ( j ) 
if (aspec (j ) 
if (aspec (j ) 
if (aspec (j ) 


•  eq. 

'Ag') 

.eq. 

'Al' ) 

.  eq. 

'B'  ) 

.  eq. 

'Ba' ) 

.eq. 

'Be' ) 

.eq. 

'Bi' ) 

.eq. 

'C'  ) 

.eq. 

'Ce' ) 

.eq. 

'Cr' ) 

•  eq. 

'Cs') 

.  eq. 

'Cu' ) 

.  eq. 

'Fe' ) 

.eq. 

'H'  ) 

.eq. 

'Hg' ) 

.eq. 

'K'  ) 

.eq. 

'Li'  ) 

.eq. 

'Mg'  ) 

.eq. 

'Mn' ) 

.eq. 

'Mo' ) 

.eq. 

'Na' ) 

.eq. 

'Ni' ) 

.eq. 

'P'  ) 

.eq. 

'Pb' ) 

.eq. 

'Rb'  ) 

.eq. 

'S'  ) 

.eq. 

'Si' ) 

.eq. 

'Sn'  ) 

.eq. 

'Th'  ) 

.eq. 

'Ti'  ) 

.eq. 

'U'  ) 

.eq. 

'V'  ) 

.eq. 

'W'  ) 

•  eq. 

'Zn' ) 

ostoic=ostoic+alp ( j ) 
ostoic=ostoic+alp ( j ) *3 . 0/2 . 0 
ostoic=ostoic+alp ( j ) *3 . 0/2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) 
ostoic=ostoic+alp ( j ) *3 . 0/2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *3 . 0 
ostoic=ostoic+alp ( j ) *3 . 0/2 . 0 
osCoic=ostoic+alp ( j ) 
ostoic=ostoic+alp ( j ) *3 . 0/2 . 0 
ostoic=ostoic+alp ( j ) /2 . 0 
ostoic=ostoic+alp ( j ) 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) 
ostoic=ostoic+alp( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *3 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) 
ostoic=ostoic+alp ( j ) *5 . 0/2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp ( j ) *2 . 0 
ostoic=ostoic+alp( j ) *5. 0/2.0 
ostoic=ostoic+alp ( j ) *3 . 0 
ostoic-ostoic+alp (j ) 


if(aspec(j)  .eq.  'Zr')  ostoic=ostoic+alp ( j ) *2 . 0 
if{aspec(j)  .eq.  ' 0' ) oin=alp { j ) 

10  continue 

equiv=ostoic/oin 

return 
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